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Abstract 

0^ ' We discuss solution concepts for linear hyperbolic equations with coefficients of regularity 

below Lipschitz continuity. Thereby our focus is on theories which are based either on a 
generalization of the method of characteristics or on refined techniques concerning energy 

i^h ' estimates. We provide a series of examples both as simple illustrations of the notions and 

, conditions involved but also to show logical independence among the concepts. 

AMS 2000 subject classification: 35D05,35D10,46F10,46F30. 

Introduction 

\^ ', According to Hurd and Sattinger in 23J the issue of a systematic investigation of hyperbolic partial 

differential equations with discontinuous coefficients as a research topic has been raised by Gelfand 
in 1959. Here, we attempt a comparative study of some of the theories on that subject which have 
been put forward since. More precisely, we focus on techniques and concepts that build either on 
the geometric picture of propagation along characteristics or on the functional analytic aspects of 
■ energy estimates. 



In order to produce a set-up which makes the various methods comparable at all, we had to 
stay with the special situation of a scalar partial differential equation with real coefficients. As 
a consequence, for example, we do not give full justice to theories whose strengths lie in the 
application to systems rather than to a single equation. A further limitation in our choices comes 
from the restriction to concepts, hypotheses and mathematical structures which (we were able to) 
directly relate to distribution theoretic or measure theoretic notions. 

To illustrate the basic problem in a simplified lower dimensional situation for a linear conservation 
law, we consider the following formal differential equation for a density function (or distribution, 
or generalized function) u depending on time t and spatial position x 

d t u(t, x) + d x (a(t, x)u(t, x)) = 0. 

Here, a is supposed to be a real function (or distribution, or generalized function) and the deriva- 
tives shall be interpreted in the distributional or weak sense. This requires either to clarify the 
meaning of the product a ■ u or to avoid the strict meaning of "being a solution" . 
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An enormous progress has been made in research on nonlinear conservation laws (cf., e.g. \17\ 12] 
and references therein) of the form 

d t u(t, x) + d x {g(u(t, x))) = 0, 

where g is a (sufficiently) smooth function and u is such that g(u) can be defined in a suitable 
Banach space of distributions. Note however, that this equation does not include linear operators 
of the form described above as long as the nonlinearity g does not include additional dependence 
on (t,x) as independent variables (i.e., is not of the more general form g(t,x,u(t,x))). Therefore 
the theories for linear equations described in the present paper are typically not mere corollaries 
of the nonlinear theories. Essentially for the same reason we have also not included methods based 
on Young measures (cf. [T71 Chapter V]). 

Further omissions in our current paper concern hyperbolic equations of second order. For advanced 
theories on these we refer to the energy method developed by Colombini-Lerner in [7] . An overview 
and illustration of non-solvability or non-uniqueness effects with wave equations and remedies using 
Gevrey classes can be found in [ST] , 

Of course, also the case of first-order equations formally "of principal type" with non-smooth 
complex coefficients is of great interest. It seems that the borderline between solvability and non- 
solvability is essentially around Lipschitz continuity of the coefficients (cf. [231 ED 122] ) ■ Moreover, 
the question of uniqueness of solutions in the first-order case has been addressed at impressive 
depth in 0. 

Our descriptive tour with examples consists of two parts: Section 1 describes concepts and theories 
extending the classical method of characteristics, while Section 2 is devoted to theories built on 
energy estimates. All but two of the theories or results (namely, in Subsections 1.3 and 2.3.2) 
we discuss and summarize are not ours. However, we have put some effort into unifying the 
language and the set-up, took care to find as simple as possible examples which are still capable 
of distinguishing certain features, and have occasionally streamlined or refined the original or 
well-known paths in certain details. 

In more detail, Subsection 1.1 starts with Caratheodory's theory of generalized solutions to first- 
order systems of (nonlinear) ordinary differential equations and adds a more distribution theoretic 
view to it. In Subsection 1.2 we present the generalization in terms of Filippov flows and the 
application to transport equations according to Poupaud-Rascle. Subsection 1.3 provides a further 
generalization of the characteristic flow as Colombeau generalized map with nice compatibility 
properties when compared to the Filippov flow. In Subsection 1.4 we highlight some aspects 
or examples of semigroups of operators on Banach spaces stemming from underlying generalized 
characteristic flows on the space-time domain. We also describe a slightly exotic concept involving 
the measure theoretic adjustment of coefficients to prescribed characteristics for (l + l)-dimensional 
equations according to Bouchut- James in Subsection 1.5. 

Subsection 2.1 presents a derivation of energy estimates under very low regularity assumptions on 
the coefficients and also discusses at some length the functional analytic machinery to produce a 
solution and a related weak solution concept for the Cauchy problem. Subsection 2.2 then compares 
those three theories, namely by Hurd-Sattinger, Di Perna-Lions, and Lafon-Oberguggenberger, 
which are based on regularization techniques combined with energy estimates. Finally, Subsection 
2.3 briefly describes two related results obtained by paradifferential calculus, the first concerning 
energy estimates and the solution of the Cauchy problem for a restricted class of operators, the 
second is a method to reduce equations to equivalent ones with improved regularity of the source 
term. 

As it turns out in summary, none of the solution concepts for the hyperbolic partial differential 
equation is contained in any of the others in a strict logical sense. However, there is one feature 
of the Colombeau theoretic approach: it is always possible to model the coefficients and initial 
data considered in any of the other theories (by suitable convolution regularization) in such a way 
that the corresponding Cauchy problem becomes uniquely solvable in Colombeau's generalized 



2 



function algebra. In many cases the Colombeau generalized solution can be shown to have the 
appropriate distributional aspect in the sense of heuristically reasonable solution candidates. 

0.1 Basic notation and spaces of functions, distributions, and general- 
ized functions 

Let fl denote an open subset of M. n . We use the notation K <s f2, if K is a compact subset of f2. 
The letter T will always be used for real number such that T > 0. We often write Ht to mean 
]0, T[xW l with closure Th = [0, T] x R n . 

The space C°°(f2) consists of smooth functions on fl all whose derivatives have continuous ex- 
tensions to f2. For any s £ M and 1 < p < oo we have the Sobolev space W s,p (M. n ) (such that 
W 0,p = L p ), in particular H s (M. n ) = W s,2 (W n ). Our notation for i? s -norms and inner products 
will be ||.|| and (., .) s , in particular, this reads ||.|| and (., .)o for the standard L 2 notions. 

We will also make use of the variants of Sobolev and L p spaces of functions on an interval J C 
R with values in a Banach space E, for which we will employ a notation as in L (J;E), for 
example. (For a compact treatment of the basic constructions we refer to [3H1 Sections 24 and 39].) 
Furthermore, as usually the subscript 'loc' with such spaces will mean that upon multiplication by 
a smooth cutoff we have elements in the standard space. We occasionally write AC( J; E) instead 
of W 1 ' c (J; E) to emphasize the property of absolute continuity. 

The subspace of Distributions of order k on £1 (k £ N, k > 0) will be denoted by X> /fe (f2). We 
identify D (O) with the space of complex Radon measures /i on O, i.e., p = v+ — v- + i(r)+ — ??_), 
where v± and r\± are positive Radon measures on f2, i.e., locally finite (regular) Borel measures. 

As an alternative regularity scale with real parameter s we will often refer to the Holder- Zygmund 
classes C*(IR n ) (cf. [TTj Section 8.6]). In case < s < 1 the corresponding space comprises the 
continuous bounded functions u such that there is C > with the property that for all x ^ y in 
K™ we have 

\u(x) -u(y)\ < c 
\x-y\ s 

Special types of distributions on M. will be used in several of our examples to follow: the Heavisidc 
function will be understood to be the L°°(M) class of the function defined almost everywhere by 
H(x) = (x < 0), H(x) = 1 (x > 0), and will again be denoted by H\ the signum function is 
sign(x) = H(x) — H(— x) ; furthermore, x + denotes the continuous function with values x+ = 
(x < 0), x + — x (x > 0), a;_ = x + — x; S denotes the Dirac (point) measure at (in any dimension). 

Model product of distributions: A whole hierarchy of coherent distributional products has 
been discussed in 28, Chapter II], each of these products yielding the classical pointwise multipli- 
cation when both factors are smooth functions. The most general level of this hierarchy is that of 
the so-called model product of distributions u and v, denoted by [u • v] if it exists. 

We first regularize both factors by convolution with a model delta net (p £ ) e> Q, where p e (x) = 
p(x/e)/e n with p 6 P(M") such that J p{x)dx = 1. Then the product of the corresponding 
smooth regularizations may or may not converge in T>' . If it does, the model product is defined 

by 

[u ■ v] = lixaOu * p e )(v * p e ). 

£—►0 

In this case, it can be shown that the limit is independent of the choice of p. For example, we 
have [H ■ S] = 8/2 and [5 ■ 5] does not exist. 

Colombeau generalized functions: Our standard references for the foundations and some 
applications of Colombeau's nonlinear theory of generalized functions are [H El [23 H2] • We will 
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employ the so-called special variant of Colombeau algebras, denoted by Q s in [H], although here 
we shall simply use the letter Q instead. 

Let us briefly recall the basic constructions and properties. Colombeau generalized functions 
on Vt are defined as equivalence classes u — \(u £ ) £ ] of nets of smooth functions u e G C°°(f2) 
(regularizations) subjected to asymptotic norm conditions with respect to e G (0, 1] for their 
derivatives on compact sets: in more detail, we have 

• moderate nets £m(^): (u £ ) £ G C°°(0)(°' 1 l such that for all K <e fl and a G N n there exists pel 
such that 

(0.1) \\d a u e \\ Lx{K) = 0(e-n (e-0); 

• negligible nets Af(il): (u £ ) e G £m(^) such that for all K <e fi and for all q G R an estimate 
KIIl~(*) =0(s q ) (e-f 0) holds; 

• £m(Q) is a differential algebra with operations defined at fixed e, Af(il) is an ideal, and Q(fl) := 
£M(^)/A/"(r2) is the (special) Colombeau algebra; 

• there are embeddings, C°°(Q) G(Q) as a subalgebra and V(Q) <^-» Q(Q.) as a linear subspace, 
commuting with partial derivatives; 

• Q —> Q(Q) is a fine sheaf and G c (fl) denotes the subalgebra of elements with compact support; by 
a cut-off in a neighborhood of the support one can always obtain representing nets with supports 
contained in a joint compact set; 

• in much the same way, one defines the Colombeau algebra Q (f2) on the closure of the open set f2 
using representatives which are moderate nets in C°°(f2) (estimates being carried out on compact 
subsets of fi); 

• two Colombeau functions u = [{u s ) e ] and v = [(v £ ) £ ] are said to be associated, we write u ~ v, 
if u e — v £ ^0 in V as e — > 0; furthermore, we call u associated to the distribution id 6 D', if 
u £ — > w in V as e — > 0; w is then called the distributional shadow of u and we also write www; 

• assume that f2 is of the form Q, = ]T\, Ti\ xf2', where fi' C R ra open and — oo < T% < T-i < oo; then 
we may define the restriction of u — [(u s ) s ] £ Q(fl) to the hyperplane {^o} X Fl' (T\ < to < I2) to 
be the element u \t=t defined by the representative (u e (to, .)) e ; similarly, we may define the 
restriction of u € G{[T\, T2] x f2') to t = to f° r 7i < to < T2 and obtain an element u \t=t G{Q')- 

• the set R of Colombeau generalized real numbers is defined in a similar way via equivalence classes 
r = [(r e )e\ °f ne t s °f re& l numbers r £ 6 R subjected to moderateness conditions |r e | = 0(e~ p ) 
(e — » 0, for some p) modulo negligible nets satisfying |r e | = 0(e 9 ) (e — > 0, for all (7); if A C R 
we denote by A the set of all generalized numbers having representatives contained in A (for all 
e S]0, 1]). Similarly, if B C R n we construct B C R n from classes of nets (x £ ) £ with x £ € B for all 

e; 

• a Colombeau generalized function m = [(u e ) e ] G ^ is said to be c-bounded (compactly 
bounded), if for all K\ <g £1 there is K 2 <s R d and £o > such that u e {K{) C holds for all 
e > e . 

1 Solution concepts based on the characteristic flow 

In this section we introduce solution concepts for first order partial differential equations, which 
are based on solving the system of ordinary differential equations for the characteristics and using 
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the resulting characteristic flow to define a solution. To illustrate the basic notions we consider 
the following special case of the Cauchy problem in conservative form 

n 

Lu := d t u + 9 *k 0*k(*. x)u) = 0, u(0) = u a G P'(M"), 
fc=i 

where the coefficients au are real- valued bounded smooth functions. The associated system of 
ordinary differential equations for the characteristic curves reads 

ih{s) = a fc (s,£(s)), £k(t) = Xk (k=l,...,n). 

We use the notation £(s;t,x) — (£i(s;t,x), ...,£ n (s;t,£c)), where the variables after the semicolon 
indicate the initial conditions x = (x\, ...,x n ) at t. We define the smooth characteristic forward 
flow 

X:[0,T]xR"^r. (s,x) ^>£(s;0,x) 
Note that \ satisfies the relation (d x denoting the Jacobian with respect to the x variables) 

d tX (t, x) = d xX (t, x) ■ a(t, x) V(t, x) G [0, r]xi", 

which follows upon differentiation of the characteristic differential equations and the initial data 
with respect to t and Xk (k — 1, . . . , n). Using this relation a straightforward calculation shows 
that the distributional solution u G C°°([0, T\; V'(M. n )) to 

Lu = 0, u(0) = u G V'(R n ) 

is given by 

(u(t),ij)~(uo,i>(x(t,.))) VfeP(M"),0<(<T. 

If there is a further zero order term b ■ u in the differential operator L, then the above solution 
formula is modified by an additional factor involving b and x accordingly. 

In a physical interpretation the characteristic curves correspond to the trajectories of point par- 
ticles. This provides an idea for introducing a generalized solution concept when the partial 
differential operator has non-smooth coefficients: As long as a continuous flow can be defined, the 
right-hand side in the above definition of u is still meaningful when we assume uq G T> (R n ). The 
distribution u defined in such a way belongs to AC([Q, T};T> °(R n )) and will be called a measure 
solution. 

This approach is not limited to classical solutions of the characteristic system of ordinary differen- 
tial equations, but can be extended to more general solution concepts in ODE theory (for example, 
solutions in the sense of Filippov) . Although such a generalized solution will lose the property of 
solving the partial differential equation in a distributional sense it is a useful generalization with 
regard to the physical picture. 

1.1 Caratheodory theory 

Let T > and Qt =]0, T[xR n . Classical Caratheodory theory (cf. [H] [Chapter 1]) requires the 
coefficient a — (ai, a n ) to satisfy 

(i) a(t, x) is continuous in x for almost all t G [0, T], 

(ii) a(t, x) is measurable in t for all fixed x G R n and 

(iii) sup^gjjn |ct(£, £c) | < P(t) almost everywhere for some positive function (3 G L 1 ([0,T]). 
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Then the existence of an absolutely continuous characteristic curve £ = (£1, ...,£„), which fulfills 
the ODE almost everywhere, is guaranteed. Note that the first two Caratheodory conditions 
ensure Lebesgue measurability of the composition s i— » a(s, f(s)) for all / E AC([0, T]) n , while the 
third condition is crucial in the existence proof. 

A sufficient condition for forward uniqueness of the characteristic system is the existence of a 
positive a £ L 1 ([0,T]), such that ((., .) denoting the standard inner product on R") 

(a(t,x) - a(t,y),x- y) < a(t)\x - y\ 2 

for almost all (t, x), (t, y) E £It (cf. [H Theorem 3.2.2]). As well-known from classical ODE theory, 
forward uniqueness of the characteristic curves yields a continuous forward flow 

X ■ [0, T] x E™ — > E™. [s,x) ^£(s;0,x) 

It is a proper map and for fixed time x(t> ■) is onto. For the sake of simplicity we assume a E 
C([0,T] x E")" and b E C([0,T] x E n ). 

Let 

i 



h b (t,x) :— exp | - / /<( r. \ ( r. r ) ) 
then u £T) (Oy) defined by 

T 

(1.2) (u,v)v'{Q T ) ■= J ( u o,tp(t,x(t,-))h b (t,-)) v ,„ (Rn) dt 



(note that it can be regarded as element in AC([0,T];T> °(E n )), so the restriction u(0) is well- 
defined and equal to uo E 2? °(E™)) solves the initial value problem 

n 

Lu := dtu + <3 X11 , (<2fc ■ it) + 6it = 0, it(0) = «o 
fc=i 

on , where a k • it and 6 • u denotes the distributional product defined by 
• : C{fl T ) x V°(n T ) -> X>'°(ft T ) 

Applying L on it we obtain 

n 

{Lu, cp) v > (0t) = (u, a A h <p + btp) v > (0t) 

(u , {-d t f - akd Xk ^ + btp)(t, x(t, -))h b (t, 0)z>'o(H»)dt. 
fc=i 

Set 4>(t,x) := ip(t,x(t,x)) and ijj(t,x) := 4>(t,x) • h b (t,x), then we have 

dt<t>(t,x) = d t <p(t,x(t,x)) = (d t ip + ^2a k (t, x)d Xk (p)(t,x(t, x)), 



k=l 
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and 



d t ip{t, x) = d t c/>(t, x)h b (t, x) + <j>(t, x)d t h b {t 1 x) 

n 

= (d t <p + ^2 ak ( t ' x ) d x k (fi)(t, x(t, x)) • h b (t, x) ~ <p(t, x(t, x))b(t, x(t, x))hb(t, x) 

k=l 

n 

= (d t (fi + ^a k (t,x)d Xk ip - bip)(t,x(t,x)) ■ h b (t,x), 

k=l 

thus 

T T 

{Lu,Lp) v , (nT) = - J {u ,d t ip(t,-)) v , 0(]&n) dt = - J d t (uo,^(t,-)) V 'o {Rn) dt = 0. 



for all (p £ T>(Ht)- The initial condition u(0) = uq is satisfied, since x(0, x) — x, thus h b (0, x) = 1. 
Remark 1.1. In this sense, we can obtain a distributional solution for the Cauchy problem 

n 

Pv := d t v + 2J a k d Xk v + cv = 0, v(0) = v , 
fc=i 

whenever a e C([0,T] x W l ) n and c e 2?'([0,T] x R"), such that -div(a) + c S C([0,T] x E n ) 
and vo G V °(K n ). We simply set b := — div(a) + c and construct the solution as above. In other 
words, such a solution solves the equation in a generalized sense, relying on the definition of the 
action of Q :— X)fe=i a kdk + c on a distribution of order by 

n 

(Qv, p) v > {nT) := -{v, akd Xk tp) v > 0{nT) - (v, (-div(a) + c)(p) v > 0{nT y 
k—l 

In case where div(a) and c are both continuous, we can define the operator Q classically by using 
the product ■ : D'°(fi T ) x C(0 T ) -> T>'°(Q T ) as above. 



1.2 Filippov generalized characteristic flow 

As we have seen in the previous subsection, forward unique characteristics give rise to a contin- 
uous forward flow. But in order to solve the characteristic differential equation in the sense of 
Caratheodory, we needed continuity of the coefficient a in the space variables for almost all t. In 
case of more general coefficients a £ Lj oc (M, L 00 (R™))™ we can employ the notion of Filippov char- 
acteristics, which replaces the ordinary system of differential equations by a system of differential 
inclusions (cf. [H]). The generalized solutions are still absolutely continuous functions. Again, 
the forward-uniqueness condition on the coefficient a 

(1.3) (a(s,x) - a(s,y),x- y) < a(s)\x - y\ 2 

almost everywhere yields unique solutions in the Filippov generalized sense. The generated Filip- 
pov flow is again continuous and will enable us to define measure- valued solutions of the PDE (cf. 
[3TJ]). as before. 

In the Filippov solution concept the coefficient is replaced by a set- valued function (t, x) — > A t ^ x Q 
1™. It has to have some basic properties which imply the solvability of the resulting system of 
differential inclusions 

£f(s) G ^4 s ,£ f ( s ), a.e., fp(i) = x, 
with £p G AC([0, oo[)™. These basic conditions are 
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(i) At jX is non-empty, closed, bounded and convex for all x E M. n and almost all t E [0, T], 

(ii) {t E [0, T] | sup a£At (a, w) < p} is Lebesgue measurable for all x E K" , w E R™ , p E R, 

(iii) for almost all t € [0, T], the set (J {x} x A ti!t is a compact subset of E™ x M" for X CC R'\ 
and 

(iv) there exist a positive function (3 E i 1 ([0, T]) such that sup agj4 |a| < for almost all 
t E [0,T] and all x E R™. 

There are several ways to obtain such a set- valued function A from a coefficient a E L 1 ( [0 , T] ; L °° (R n ) ) " , 
such that the classical theory is extended in a compatible way. Thus the corresponding set-valued 
function A should fulfill A to . Xa '■= {a(to, %o)} whenever a is continuous at (to,xo) E [0, oo[xR™. 

One way obtaining a set-valued function corresponding to a a E L 1 ([0, T]; L°°(R n )) n is by means 
of the essential convex hull ech(a). It is defined at (t, x) E [0, T] x R n by 

(ech(o))t, x := P| P| ch{a{t,B s (x)/N)) 

A(iV) = 

where ch(M) denotes the convex hull of a set M C R™ and A is the Lebesgue measure on R n . 

Another way is to use a mollifier p E 5(R") with J p(x)dx = 1, put p e (x) — e~ n p(e~ 1 x) and 
A £ a * p e | [0 ,t]xK" where a E L°°(R" +1 )" is the extension of a function a E L°°([0, T] x R n ) n 
by zero. Then the concept of a generalized graph Ca as defined in [13J yields a set-valued function 
satisfying the above basic properties. 



1.2.1 Measure solutions according to Poupaud-Rascle 



Let floo := ]0, oo[xl". We assume a E L 1 1 oc (IR + ; L°°(W l )) n to be a coefficient satisfying the forward 
uniqueness criterion p.3|) . Let Lu := dtu + Y^i=i 9 Xi (ciiu) an d £f be the unique solution to 

(1.4) £ F (s) £ ech(a) S)?J?(s ), = Z. 
The map 

X f ■ R+ x R" -> R™, (i,ar) i-> £ F (t;0,x) 
is the continuous Filippov (forward) flow. 

Definition 1.2 (Solution concept according to Poupaud-Rascle). Let uq E Al(,(M) n &e a bounded 
Borel measure, then the image measure at t E [0, oo[ is 

(1.5) «(*)(£):= / 1b(xf(M))^oW, 

w/iere B C R" is some Borel set. The map u: [0,oo[-> M b (W 1 )) belongs to C([0, oo[; M b (M. n )) 
and is called a measure solution in the sense of Poupaud-Rascle of the initial value problem 

n 

Lu := d t u + }^ d Xk (ofe ■u) = 0, u(0) = w . 
iVo£e £/ia£ u defines a distribution of order m P'^oo) &y 

(«, v)-D'(n oc ) := y ("o, v(t, XF(t, x))) V 'oQgp-) dt, V<p E P'(fioo)- 
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The solution concept of Poupaud-Rascle does not directly solve the partial differential equation in 
a distributional sense, but it still reflects the physical picture of a "transport process" as imposed 
by the properties of the Filippov characteristics. Nevertheless, in the cited paper of Poupaud- 
Rascle(j3U]) the authors present an a posteriori definition of the particular product a ■ u, which 
restore the validity of the PDE in a somewhat artificial way. We investigate this in the sequel in 
some detail. 

Definition 1.3 (A posteriori definition of a distributional product in the sense of Poupaud-Rascle). 

Let u G T>'(W l ) be a distribution of order and a G L 1 1 oc ([0, oof, L°°(R n )) , satisfying the forward 
uniqueness condition 11. 3\) . such that there exists a continuous Filippov flow \F- Furthermore we 
assume that u is a generalized solution of the initial value problem as defined in 11. 5\) . Then we 
define the product a • u — (a/. • u)k in T>'(]0, oc[xR™)™ by 



(a«ii,^)x)'(noo) : = ( u o, j dtXF(t,x)(p(t,x(t,x))dt) T>l0{Rn) , (peV(n x ). 

o 

Remark 1.4. Note that the product a-u is defined only for distributions u that are a generalized 
solutions (according to Poupaud-Rascle) of the initial value problem (II. 6|) with the coefficient a. 
The domain of the product map (a, u) ^ a»u, as subspace of T> °(R") x T> °(R n ) has a complicated 
structure: Just note that the property to generate a continuous characteristic Filippov flow \f 
is not conserved when the sign of the coefficient a changes, as we have seen for the coefficient 
a{x) — sign(x). 

Example 1.5. Consider problem (|f .6p with the coefficient a(x) := — sign(x) subject to the initial 
condition uq = 1. Then the continuous Filippov flow is given by 

XF(t, x) = -{t + x)-H(-x) + (x- t)+H{x). 

We have XfQ, 0) = t+ - (-*)_ = and 

dtXF(t, x) = -H(-t - x)H(-x) - H(x - t)H(x) for almost all t e [0, oo[. 

The generalized solution u is defined by (u, ip) := (uq, <p(t, x)) dt, where (f>(t, x) := tp(t, \{t, x)). 
We have that 

( ip(t,x + t) x < 0,0 < t < -x 
4>(t,x) := I tp(t,Q) t > \x\ 

[ (p(t,x- t) x > 0,0 < t < x, 

thus 



u > V)i5'o(o 00 ) : = / {uQ,4>{t,x))dt = I I 4>{t,x)dx dt 



-oo 

oo oo — t oo oo 

2 J (p(t,0)tdt + J J tp(t,x + t)dx dt + J J (p(t,x-t)dx dt 

-oo t 

oo oo oc 



= 2 J ip{t,0)tdt + J J <p(t,z)dz dt + J J tp(t,z)dz di= <l + 2*<y,<p(t,-))z>'o( K » ) 

-oo 

This generalized solution gives rise to the following product 

((-sign(a;)) • (l + 2t6(x)),<p) T>l (n oo) := (1, / d t ^F(t,x)ifi(t,^ F {t,x))dt) v ' 0(Mn) 

o 
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in D'(Ooo). Evaluating the right-hand side we obtain 

oo oc / oo 

(l,jd tX F(t,x)<p(t,XF(t,x))dt) = J I- J H(-x)H(-x-t)<p(t,x + t) 



OO 

J H(x)H(x 



t)ip(t, x — t)dt dx. 



Since H(—x)H(—x — t) = H(—x — t) and H{x)H(x — t)= H(x — t) for t > the latter gives upon 
substitution 

oo oo oo 

(1, J dtXF(t,x)tp(t,XF(t,x))dt) = - J sign(z) y ip(t,z) dtdz, 

-oo 

hence 

(— sign(x)) • (1 + 2tS(x)) = — sign(x). 

However, we cannot define the product if — sign(x) is replaced by +sign(x), since the Filippov 
characteristics £f(£;0, x) are no longer forward unique and thus do not generate a continuous 
Filippov flow xf- 

Example 1.6. We consider the same coefficient a(x) :— — sign(x) as before, but now we set 
u := S. We obtain the generalized solution 

oo 

(".vVcnoo) := { l ®s,v{t,XF{t,x))) v , (noo) = J ip{t,XF{t,Q))dt 



This enables us to calculate the product 

oo 

((-sign(x)) •5{x),<p) = -{6, J dtXF(t,x)<p(t,£ F (t,x))dt). 

o 

Putting ij){x) = J °° dt£, F {t, x)ip(t, £ F (t, x))dt and observe that 

oo —X 

V>(x) := J dtXF(t,x)(p(t,XF(t,x))dt = J <p(t,x + t)dt, if x < 
o o 

and 

oo X 

ip(x) — J dtXF(t,x)tp(t,XF(t,x))dt = - J ip(t,x — t)dt, if x > 0. 



At x = we obtain i[)(Q) = lim^^o- "4>(x) = linia;^o + ^P( x ) = 0, so it follows that (—sign) • S = 0. 
Example 1.7. Let a(t,x) := 2H(—x), so that the Filippov flow is given by 

X F (t, x) = -(x + 2t)-H{-x) + xH(x). 
We have xf(*, 0) = -2i_ = and 

d t x F (t, x) := 2H(-x - 2t)H(-x). 
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Hence dtXF^i 0) = for almost all t <G [0, oo[. If uo = 1 the generalized solution is 

oo oo 

( u ^)v'"(n x ) '■= J J <P(t,x)dxdt, 

-oo 

where <p(t,x) — (p(t,XF(t,x)). Since 



ac) |{x<-2t} = x + 2t) 

<f>(t,X) \{~2t<x<0} = V(t,0) 

<f>(t,x) \{ 0<x } = (p(t,x), 



we obtain 



(u,<p) T) >0(n oo )= J I / ¥>(t,x + 2t)dx + 2tip(t,0) + J ip(t,x) dx ] dt = <l+t<5, ^(i, -))x>'o(r»)! 

\-oo 

hence u = 1 + 2i<5(x). Again we determine the product (2iJ(— a;)) • (1 + 2t5(x)) by 

oo oo 

(2H(-x) • (1 + 2*5(x)), = 2 / / H(-x)H(-x - 2t)<p(t, x + 2t) di dx 



-oo 

oo oo oo oo 

= 2 y J H(-x-2t)<f(t,x + 2t)dtdx = 2 J J H(-z)ip(t, z) dz dt = (1 ®2H(—),(p} v ,( Qoo) . 

— oo — oo 

We obtain (2H(— x)) • (1 + 2t5{x)) = 2H(—x). Observe that together with the result in Example 
(fL5| (-sign(x)) • (1 + 2i(5(x)) = (2H(-x) - 1) • (1 + 2i<5(x)) we can conclude that either (-1) • 
(1 + 2t5(x)) is not defined or the product • is not distributive. In fact , it is not difficult to see 
that (—1) • (1 + 2tS(x)) cannot be defined in this way, neither can 1 • (1 + 2tS(x)). 

Example 1.8 (generalization of Example II. 5ft . Let c\ > C2 be two constants, and a € [ci,ca]. 
Consider the a(t, x) := ciH(at — x) + C2H(x — at). We set t\{x) := c ~* a if x < and t2{x) :— a ^ C2 
for x > 0. The unique Filippov flow is given by 

c\t + x x < 0, t < ti (x) 
at x < 0, t > ti(x) 

XF(t,x) = < oct x = 0, 

C2t + X X > 0, t < t2{x) 

at x > 0, t > t 2 (x) 



The generalized solution of the initial value problem Lu := dtu+d x (a-u) = 0, u(0) = uq 6 Lj oc (l 
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according to Poupaud-Rascle is given by 

T ti(x) 

(u, <^)-p'o (0t) = J {u , <p(t, xp(t, -)))x)'o(R) dt = J J u o{x)ifi(t, ctt + x) dt dx 

— oo 

T 00*2(2;) 00 T 

+ J J uo(x)ip(t, at) dt dx + J J uo(x)tp(t, c$t + x) dt dx + J J uo(x)<p(t, at) dt dx 

-DOtl(z) 00 t 2 (x) 

T -t(ci-a) T 00 

= J J uo(x)tp(t, c\t + x) dx dt + J J uo(x)ip(t, C2t + x) dx dt 

-00 -t(c 2 -a) 

T / t(a-c 2 ) 
+ J I J Uq(. <■)'!.! J Y *il.nl)<ll. 

\-t( Cl -a) 

hence 

/ t(Q-c 2 ) 

u := uq(x — c\t)H(at — x) + uq(x — c%t)H{x — at) + J uq(x) dx |j ,; - n / ). 

\-t(ci-a) 



1.3 Colombeau generalized flow 

In this subsection we consider the solvability of the ordinary differential equations for the charac- 
teristics in the setting of Colombeau generalized functions. Our main focus will be on distributional 
shadows of such generalized solutions. It will appear that under certain assumptions on the right- 
hand side, the distributional shadow exists and is absolutely continuous. We will also show a 
uniqueness result for distributional shadows. 

Theorem 1.9 (Existence). Assume A 6 Q(Qx) n with a representative (A £ ) £ , such that 
(1.6) sup |A e (t,x)| < f3(t), e e]0, 1], almost everywhere intE [0,T] 

holds, where (3 is some positive function in i 1 ([0, T}). Let (t,x) £ 12t be the initial condition. 
Then there exists a c-bounded solution £ G C/([0,T]) n to the initial value problem 

i(s) = A(s,t(s)), £®=x. 

Furthermore, there exists some (t, x) s £,c £ ^-C([0, T]) such that for any representantive (£ E ) E 
of £, (t £ ,x £ ) £ of(t,x) there exists subsequences (t £j , x e .)j,(£ £ .)j withlimj^ oc (t £j , x £j ) = (t,x) and 

t; £j 3 — ► t;c uniformly on [0, T] and £c(t) = x. 

Proof. By classical existence and uniqueness we obtain £ e for each e s]0, 1] such that 

s 

&(*) = X £ + / A £ {T,£ £ {T))dT 



holds. Condition (|1.6p yields |£ £ (s)| < |a? e | + | J* 0(t) dr\ for all s,t 6 [0,T], hence c-boundedness 
of (£ e ) e on [0,T] and furthermore moderateness of £ e (by [13] [Proposition 1.2.8]). In fact this 
existence result is quite similar to the one given in [^[Proposition 1.5.7]. 
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To prove the existence of a convergent subsequence of (£ e ) e , we may assume without loss of 
generality that lim e ^o(ie, %e) — (t,x) G fi-r- Note that the family (£ E ) e is uniformly bounded and 
equicontinuous, since 

s' 

|&(«) -&(«')!< I Jf3(r)dr\ s,s'e [0,T],e€]0,l]. 

s 

The Theorem of Arzela-Ascoli yields a subsequence (£ e ^ )j converging uniformly to some £c G 
C([0,T]). Clearly, £ c (t) = lim^oo £ Sj (t e J = lim i _ ) . 00 x Ej = x. such that 

lim sup (s) - £c(s)| = 0, 

J~»oO;5 e [0,T] 

holds. We have for all s, s' 6 [0, T], 

< -6,(3)1 + | /3(r)dT+ -£c(*')l | /3(r)dr, 

s s 

hence £c is absolutely continuous on [0, T]. □ 
1.4 Semigroups defined by characteristic flows 

Let X be a Banach space and (S t ) tg [ ,oo[ be a family of bounded operators E t on X. Consider 
the following conditions: 

(i) S = id 

(ii) S s o £ t = S s+f for all s, i G [0, oo[ and 

(iii) the orbit maps 

a u „ : [0, oo [ -> X 

t i * E t (uo) 

are continuous for every mo £ I. 

If (i) and (ii) are satisfied, then we call (St) a semi-group acting on X. If in addition property 
(iii) holds, we say (St)te[o,oo[ is a semi-group of type Cq. 

We briefly investigate how the solution concepts discussed in subsections 1.1 and 1.2 fit into the 
picture of semi-group theory when the coefficient a is time- independent. First we return to the 
classical Caratheodory case: Let a G C(R™)™ and assume that a suffices the forward uniqueness 
condition (|1.3|) . This implies that the characteristic flow \ : — > 1" is continuous and x(t, •) is 
onto R" for fixed t G [0, oo[. Furthermore we have x( s , x( r : a?)) = xi s + r J a; ) f° r a U 1 £ ^™ an d 
r, s G [0, T] with s + r G [0, T], since a is time independent. 

Consider the inital value problem Pu = <9 t + Sfc=i a kd Xk u — with initial condition u(0) = uq G 
Co(M") (i.e. vanishes at infinity). It is easy to verify that 

£ t :Co(M") C (M n ) 

M0 ^ X*M0 
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defines Co semigroup on the Banach space Co(Qoo): Note that E t is a bounded operator on Co(R n ) 
for each t £ [0, oo[, as R") = R n 5 so 

||St(uo)||oo = sup \\u (x(t,x))\\ = sup \\u (x(t,x))\\ = sup ||u (x)|| = ||u ||oo- 

a;GR™ a;GR™ 

We have that ||S t || = 1 for all t E [0, oof. Condition (i) and (ii) follow directly from the flow 
properties of x- The continuity condition (iii), which is equivalent to 

lim ||£t(uo) - «o||oo = lira sup \\u (x{t,x)) - u (x)\\ = 0, 

holds, since (x(i, x)) X £R" is an equicontinuous family and uq vanishes at infinity. 

Remark 1.10. For a coefficient a in L°°(M. n ) n we can also define a semi-group on Co(M n ) by 
S t (uo) := Uo(xF(t,x)), where xf is the generalized Filippov flow as introduced earlier. This is 
due to the fact, that the Filippov flow has almost the same properties as the Caratheodory flow. 

It seems natural to understand the solution concepts as defined by (|1.2[) and (|1.5[) as action of the 
dual semigroup (SJ?) on 

(u(t), v)i.'0( R n) = (S> , ^>Z>'0(K») = ( u 0,S t (^)) x ,' O(M „ ) . 

the Banach space of finite complex Radon measures, the dual space of Co(R") (cf. Chapter 4], 
[29} Chapter 1.10] or [371 Chapter IX. 13] for the general setting). However, in general the dual 
semi- group is not of class Co (cf. [HI Example 1.31]). This is only guaranteed if we start from a 
Co semi-group defined on a reflexive Banach space. 

Nevertheless the solution concepts in (|1.2|) and (|1.5p still yield the semi group properties (i) and 
(ii) with weak-* continuity replacing the strong continuity property (iii). 

The situation is much easier with Hilbert spaces, of course. We conclude with an example involving 
a discontinuous coefficient. 

Example 1.11. Let a £ L°°(M.) such that there exist cq,ci > such that c\ < a(x) < C2 almost 
everywhere. We want to solve the initial value problem 

Pu = d t u + a(x)d x u = 0, u(0) = u e L 2 {R) 
for u E AC([0,T};L 2 (R)) n ^([0, T]; i^O^- 
Lct A(x) = Jq a(y) _1 dy, which is Lipschitz continuous and strictly increasing (thus globally 
invertible) and observe that x(i, x) = A~ x (t + A(x)) defines the (forward)characteristic flow that 
solves 

t 

X(t,x)=x + J a(x(r,x))dT. 


Let Q := —a(x)d x with domain D(Q) := H 1 (M.). The resolvent of Q for Re(/i) > is obtained 
from the equation 

(-Q + H)v = f, feL 2 (R). 

Upon division by a we deduce 

(1.7) d x v + ^v= f -. 

a a 

Let us first consider uniqueness: Let w E H 1 ^) satisfy 

(1.8) d x w + -w = 0. 

a 
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Since w is absolutely continuous we have 

w(x) = C exp I — 2Re(^) f — !— - dz 
\ J a[z) 

\ — oo 

for some constant C. But w 6 L 2 (R) if and only if C = 0, thus w = 0. 
Existence: One easily verifies that 



v(x) = (R(n)f)(x) 



exp — u / a(z) dz — — 



exp(-MA(x)-A(y)))44^ 



is a solution of p.7[) in AC(R). Upon substitution y i— ► z = A(a;) — A(y) in the right-most integral 
we obtain 

oo 

(1.9) (R{n)f)(x)= [exp(-»z)f(x(-z,x))dz, 

o 

which is the Laplace transform of f(x(—-,x)). 

We denote the kernel of the integral operator by 

M(x,y) := H(x~y)exp(-fi(A(x)-A(y))a{y)- 1 . 

We briefly sketch the derivation of L 2 estimates for the operator powers R(^) k for Re(/i) > 0: 
Note that R{jj) k is an iterated integral operator of the form 



R(j*rf{x) := fl(At)*- 1 / M(-,z 1 )f(z 1 )dz 1 (x) 



M(x, Zk)M(z k , 2fc_i) . . . M(z 2 , zx)f(zi) dz k -i ...dz 2 dz\ 



To simplify notation let z = (zi, . . . , Zk), d k z — dz\...dzk, h(z) :— exp (— A*X)/=i z l)i an d 
g(x, z) :— f(x(— X)f=i z h x ))- Using the flow property of x we obtain that 

R(ti) k f(x)= J h{z)g{z,x)d k z 

[0,oo[ fc 

holds, hence by by the integral Minkowski inequality 

2 x 1/2 



(1.10) \\R^) k f\W < 



\h(z)\\g(z,x)\ d k z dx 



\R \JO,oo[ fc 



< 



1/2 

2 J™ I jfe . 



h{z)\ / \g{z,x)\ z dx d k z. 



Since 



k 

J \g(z,x)\ 2 dx = J \f( X (-J2z l ,x))\ 2 dx = J \f(y)\ 

TO TIE l — l wa 



dy<-\\f\\l^ 
Co 
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we conclude 



V c o 



| i2 - / \h(z)\d«Z=-±— ] ;-\\f\\ L2 . 

Re u 

]-oo,0] fc 



The Hille-Yosida theorem Theorem 5.2]) yields that Q generates the Co semigroup 

S t : L 2 (R) -» L 2 {R) 

uo i-> u (x(-t,x)). 

The resolvent operator fj, i— » R(fi) (defined for ReQu) > 0) is the Laplace transform of the semigroup 
t — ► S t as indicated in (11. 9j) . 

Remark 1.12. Since L 2 (IR n ) is reflexive the dual semigroup is Cq as well and has as its generator 
the adjoint operator Q*. 

Remark 1.13. If we assume additional regularity on the coefficient, e.g. a 6 C%(M) with a > 0, 
in Example 11.111 then we obtain a Co semigroup (S t ) acting on the Hilbert space H S (M) with 
< s < a. We may then use the fact that the (square of the) Sobolev norm \v\ s is equivalent to 
the following expression (cf. [TBI Equation (7.9.4)]) 

\v{x)-v(y)\ 2 



\v(x)\ z dx + C s 



dx dy, 



where the constant C s depends only on the dimension n and s. From this it can be shown that 
we may have D(Q) = H S+1 (M.) as domain of Q (this also corresponds to the special case of the 
mapping properties stated in |35[ Chapter 2]). Clearly uniqueness in the characteristic equation 
(|1.8|) is still valid. A corresponding variant of the estimate (|1.10j) for the powers of the resolvent 
operator R(n) k on iJ s (R) is obtained by the following calculation (with the notation h and g as 
in Example II. lip : 



\\R(v) k f\\ 2 8 = \\RM k f\\ 2 



\(R(v) k f)(x)-(RQ,) k f)(y)\< 



\x - y\ 



l+2s 



dx dy 




\h(z w 9{x : z) ~ 9 ^: z)l dz\ dxdy 



x-y\ 



g(x,z) - g(y,z)\ 2 



\x-y 



l+2s 



1/2 

dxdy I d k z 



To carry out the x and y integrations we use the substitutions x' = %(— Ez=i Zl ' x )> v' = 
X(~ Ya=i z hV) to obtain 

l{x,z) - g(y,z)\ 2 



\x~y\ 1 + 2s 



dxdy= I I l/(x(-Sf=i^.»))-/(x(-E!Li^,i/))l 2 (fa;dy 



\f(x')-f(y')? 



R E 



lx(Ef=i^^'))-x(Eti^>y , )l 1+2s 



a (x(Ef=i^>0) 



a[x') 



dx' dy' . 



Now, by the mean value theorem we have |x(E?=i z '' x )) ~ x(E/=i z l> v)\ — ffl 2 - ~ v\ an< ^ the 

< — , thus we arrive at 



°(x(Eti^,-)) 

a(-) 



assumed bounds for a give 

\g(x, z) -g(y, z) | 2 



B R 



F -y 



l+2s 



16 



Again by Jj oo[fc \h(z)\d k z = Rc ^ )k we conclude 

3+2s , \ 3+2s 



i.e., ||i?(/i) fc /l| s < VC ° 7 - - k ||/| 



(3+2s)/2 

I 

Re(/i) A 



1.5 Measurable coefRcients with prescribed characteristics 

This subsection discusses a solution concept according to Bouchut- James ([!]), which is settled in 
one space dimension and — from the distribution theoretic point of view — can be considered as 
exotic. The basic idea is to interpret the multiplication a ■ u occurring in the partial differential 
equation as a product of a (locally finite) Borel measure u and a function a from the set B°° of 
real bounded and Borel measurable functions. 



Multiplication of Radon measures by bounded Borel functions: We may identify locally 
finite Borel measures on R with (positive) Radon-measures, that is the non-negative linear func- 
tionals on C C (R) ([14, Remark 19.49]). Moreover, the space T> °(R) is the space of complex Radon- 
measures, which allows for a decomposition of any m£P °(R) in the form u = V-i r —v^+i(r)-\-—r)-"), 
where v + , v_ , rj + , r/ + are positive Radon- measures. 

The product of a bounded Borel function a E B°° (R) with a positive Radon measure fj, is defined 
to be the measure given by 

(a //)(£?) := J lB(x)a{x) d^i(x), 

R n 

for all Borel sets B in R. Clearly, a n is again a locally finite Borel measure. 

The product employed in [5] is the extension of to £>°°(R) x T> °(R) in a bilinear way, i.e. 

o : £°°(R) x D'°(R) -> 2>'°(R) 

(a, u) i > fl + 0y + +d-0i/-- (a_ v + + a + © vJ) 

+i(a + i] + + a_ ?7_ - (a_ r; + + a+ 7?-)). 

Consider the following sequence of maps: 

C b (R) ^ £°°(R) A Lg^R) ^ X>(°>'(R) 

where l\ , t2 are the standard embeddings and A sends bounded Borel functions to the correspond- 
ing classes modulo functions vanishing almost everywhere in the sense of the Lebesgue measure. 
Although we may identify C b (R) and L™ c with subspaces of £>'°(R) this is not true of B°°(R), 
since A is not injective. Note that 12 o A o ^ is injective though. The following example illustrates 
some consequences of the non-injectivity of the map A for the properties of the product o. 

Example 1.14. Let a E R and a a (x) := 1,2; 7^0 and a Q (0) := a and u = 5 E P'°(R). Note that 
\oLi(a a ) = 1 as a distribution and the standard distributional product gives Aoi 1 (a Q )-u = 1-u = 5 
for all a E R. On the other hand a a o 5 = ad. 

The product o will be used in the solution concept for transport equations on Qt = [0, T] x R with 
coefficient a in B°°(ilT) and solution u E B°°([0, T}; T> °(R)), i.e., u is a family of distributions 
(u(t))te[o.T] such that (u(t) , cp) v ' ^ is a bounded Borel function on [0,T] for all (f E C C (R). The 
extension of the product o to this space causes no difficulty. 
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The solution concept according to Bouchut and James: A key ingredient for the solution 
concept according to Bouchut- James, is to stick to a particular representative of the coefficient (in 
the L°° sense), by prescribing the value of the coefficient a at curves of discontinuity. We refer to 
the following requirements on the coefficient a G B°°(£It) as Bouchut- James conditions: Assume 
there exists a decomposition £It = CUPU5 such that 

(i) S is a discrete subset fir, 

(ii) C is open, a is continuous on C, 

(hi) I? is a one-dimensional C^-submanifold of fix, i.e., for each (tg,xo) G T> there exists a 
neighborhood V of (to,Xo) and a C 1 parametrization of the form t i— » (t, £(t)) in 2? n V. 
Furthermore, a has limit values for each (t, x) G Z? from both sides in C \ T>. These limits 
are denoted by a+{t, x) and a_(t, 2c). 

(iv) a(t, x) G [a_(t, x), a+(t, jc)] for all (t, ac) G T>, 

(v) for any point (to,xo) £ 2? with neighborhood and local parametrization £ as in (iii), we 
have = a(t,£(t)). 

Condition (v) prescribes the values of the coefficient ait, x) on the curves of discontinuity in such a 
way that the characteristic differential equation holds. In this sense, a coefficient satisfying (i) — (v) 
is a piecewise continuous bounded function, where the (non-intersecting) curves of discontinuity 
can be parametrized as regular C 1 curves. 

The Bouchut- James solution concept interprets hyperbolic Cauchy problems in (1 + 1) dimension 
as 

(1.11) Pu := d t u + aod x u = 0, u(0) = u £ BVi oc (R) 
and 

(1.12) Lu := d t u + d x (ao u) = 0, u(0) = u G V°(R). 

Note that P (resp. L) is well-defined on the set B°° ( [0, T] ; BVi oc (M) ) (resp. B°°([0, T];V°(R))). 
Now the main results of Bouchut-James [3] are: 

Theorem 1.15. Theorem 3.4] Assume that a satisfies the Bouchut-James conditions (i)-(v). 
For any u G BV loc (R) there exists u G Lip([0,T];iJ, c (R)) n B°°([0,T]; BViocfM)) sofomy j77771J 
and swc/i iftai /or any x\ < xi we have for all t G [0, T] 

Var/(u(t, •)) < Varj(wo), 

IK«,-)IU«>(7) < IKIU~(j) 

where I :=]x\,X2[ and J :=]x\ — \\a\\ oa t,X2 + ||a||oot[. // in addition the coefficient a satisfies the 
one-sided Lipschitz condition 

(a(t, x) — a(t, y), x — y) < a(t)\x — y\ 2 for almost all (t, x), it, y) G £It, 

where a G L 1 ([0, T]) 7 then the solution u is unique. 

Theorem 1.16. [3, Theorem 3.6] Assume that a satisfies the Bouchut-James conditions (i)-(v). 
Then it follows that for any u G V°(R) there exists u G C([0,T];X»'°(R)) solving [TUty . If a 
satisfies in addition the one-sided Lipschitz condition 

(a(t, x) — a(t, y),x — y) < a(t)\x — y\ 2 for almost all (t, x), it, y) G ilr, 

where a G L 1 ([0,T]) 7 then the solution u is unique. 
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We compare the solution concept of Bouchut- James with the generalized solutions according to 
Poupaud-Rascle. 

Example 1.17. We come back to Example 1 1.8[ where ait, x) := c\H{at — x) + ciH{x — at) with 
C2 < c% and a £ [c2,ci]. Let A £ i 1 ([0,T]) such that c\ > X(t) > C2- Consider a representative a 
in £>°°(f2T) of the coefficient a given by 

!C\ x < at 
X(t) x — at . 
C2 x > at 

We investigate wether the distribution u given in Example 1 1 . 81 solves (|1.12[) in the sense of Bouchut- 
James. 

Let us consider the case uo = 1. Then we obtain the solution u = 1 + t{c\ — C2)S(x — at). Note 
that the requirement that a fulfills the Souchut- James conditions (i)-(v) forces X(t) := a for all 
t e [0,T]. Thus we have 

a o v — c\H(at — x) + a§(x — at)t(c± — C2) + CiH(x — at) 

and therefore d x (aov) = at(c\ — 02)6' (x — at) — (ci — C2)S(x — at). Since dtv — (ci — C2)5{x — at) — 
at{c\ — €2)6' (x — at) we deduce that u solves the Cauchy problem in the sense of Bouchut-James. 

Finally, we check whether the differential equation is fulfilled, if we employ the model product (cf. 
[28} Chapter 7] or the introduction) instead of o. Let [a ■ v] denote the model product. We have 
that 

[a • v] = [{c\H{at — x) + C2H(x — at)) ■ u] 

= c\H(at — x) + C2H(x — at) + c\t{ci — C2)[H(at — x) ■ 5{x — at)] 

— C2t(c\ — C2)[H{x — at) ■ d(x — at)] + c\H{at — x) + C2H(x — at) + — (ci + C2){c\ — C2)5{x — at), 
hence 

d x [a-v] = (c 2 - ci)S(x - at) - ^(ci + c 2 )(ci - c 2 )5'(x - at), 
dtv — (ci — C2)5{x — at) — at(c\ — C2)5'(x — at). 

Therefore v solves the initial value problem 

d t v + d x [a-v] =0, v(0) = l, 

if a = |(ci + c 2 ). 

Note that the coefficients H(—x) (with ci = 1, c% = 0, a = 0), —H(x) (with c\ = 0, C2 = —1, a = 
0), and — sign(cc) (when c\ = 1, ca = — l,a = 0) are included as special cases of the example 
presented here. In case the coefficient reads 

a(x) := H{-x) 

the unique solution in the sense of Bouchut-James is given by 

u = 1 + tS. 

It has been shown in (2H1 Theorem 5] that no distributional solution exists in this case when the 
model product is employed. 
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2 Solutions from energy estimates 



2.1 Direct energy estimates 

We briefly review the standard techniques of energy estimates for the initial value problem 

n 

Pu := d t u + V a,j d x u + cu = f in]0,T[xR", 

(2.13) j=i 

u(0) = ti e L 2 (R n ). 

Let q G [2,oo]. We assume that / G L 1 ([0, T]; L 2 (M™)), a = (a u ...,a n ) G ^([0, T]; ff 1 ^!"))" 
with real components, c £ i 1 ([0, T]; L 9 (R™)) and in addition 

(2.14) Idiv»-c G L 1 ([0, T]; i°°(M™)). 
2.1.1 Example derivation of an energy estimate 

We browse through the typical steps that lead to an estimate in the norm of L°°([0, T]; L 2 (R n )) 
for any 

u G AC([0,T];L 2 (K n ))nL oo ([0,T];W 1 ' J, (]R")) 
with p G [2, oo] such that | + | = ^ in terms of corresponding norms for m(0) and Pu. 
We write P = d t + Q with Q := X)fe=i t)d Xk + c(x, t) and observe that 

Pu G L 1 ([0,T};L 2 (R n )) 

holds since d t u G I^QO, T]; L 2 (K™)) and Qu G i^QO, L 2 (R™)) (the latter follows from the facts 
that d Xj u(t,.) G L 2 and L q ■ L p C 1? when 2/p+2/q = 1). Hence r i-> Re(((Pu)(r), u(r)) ) is 
defined and in L 1 ([0,T]). Furthermore, the map 1 t— > ||u(i, -)|| is continuous. 

We put 

r 

ft(r):=||^div x (a(r,-))-c(r,-)||oo and \(r) := 2 J h(s) ds > (re[0,T]). 

o 

By assumption, ft G ^([O.T]) and A £ AC([0,T]). 

The standard integration by parts argument gives the Garding-type inequality 

(2.15) ^((Qu(t),u(t))o + (u(t),Qu(t)) q ) - Re((Qu(r), «(t)) ) > -ft(r)|| U (r)|| 2 , 



and thus 



T 

y e- A ( r )Rc(((Pu)(r),u(r))o)dr 
o 

r r 

= ^ J e- x ^±-\\u{r)\\ldr + J e- x ^R C {{Qu){r),u{r))v dr 
o o 

> \e-m\\u<r)\\l Ilium 2 - J Lr) &\ e-^)|| tt(r) ||gdr. 



>- 

=0 
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Therefore 

T 

e- X ^\\u(r)\\l < \\u(0)\\l + 2je-^\\(Pu)(r)\\o\\u(r)\\ dr 

o 



< |K0)||g + 2sup (e- A(r)/2 h(r)|| ) / e-^\\(Pu)(r)\\ Q 



dr. 



o 



where we may take the supremum over r G [0, t] on the left-hand side and thus replace r by f on 
the right-hand upper bound. A simple algebraic manipulation then gives 



sup \\e~ x ^/ 2 u(r)\\o- / e- x ^/ 2 \\(Pu)(r)\\ dr < \\u(0)\\ + / e- x ^/ 2 \\(Pu)(r)\\ Q dr 
K rem I 1 \ o 

Upon removing the squares and multiplying by exp(A(t) /2) we obtain the following basic inequality. 
Energy estimate: 



(2.16) sup ||u(r)|| < exp( / h(a)da) ■ ||u(0)|| + 2 exp( / h(cr)da) ■ / \\(Pu)(r)\\ dr 

re[0,t] J J 



t I t 

= exp(/ h(a)da) I ||«(0)|| + 2 / \\(Pu)(r)\\ dr 



We recall that the exponential factor depends explicitly on the coefficients a and c via h(r) = 
\\±div x {a(r, •)) - c(r, -)||oo- 

Note that this derivation of an energy estimate relied on the Garding inequality (|2.15[) . 

Example 2.1 (Failure of the Garding-inequality (|2.15[) ). Let a e ]l/2, 1[ and define a: M — > K by 
a(x) := 1 + x°l_ when x < 1, and a(x) := 2 when x > 1. We have a G C"(K) \ Lip(M). 

Let Q: iJ x (R) -> i 2 (R) be the operator defined by (Qv){x) := a(x)v'(x) for all v G F^R). Note 
that compared to the general form of the operator Q in the derivation of the energy estimate 
above we have here c = 0, a G C°° ( [0, T] ; W^ 1 ' 2 (M) ) but diva/2 -c = a'/2 ^ L°°(M). 

Since Q is time independent, inequality (|2.15[) with some h G L 1 ([0,T]) (not necessarily of the 
form given above) would imply 

3C 6 R, Vn 6 C C °°(]R) : Re({Qv,v) ) > -~C\\vg. 

We will show that there is no constant C G R such that the latter holds. Thus (|2.15p cannot hold 
for Q (for any h G L 1 ^^})). 

Let p G C 1 (R) be symmetric, non-negative, with support in [—1,1], ||/d[|o = 1, and such that 
p'{x) < when < x < 1. We define v e (x) := e~ 1/2 p(x/e) (x G R, £ > 0). Then clearly 



21 



v E G C C °°(R) C tf^K) and ||v e || = 1 for all e > 0, but 
(Qv e ,v e )o= a(x)v' e (x)v e (x)dx 



1 oo 

v' E (x)v e (x) dx + J x a v' e {x)v e {x) dx + J v' e (x)v e (x) dx 
o 1 



=o =o 

l 

= e a ~ l J z a P '{z)p{z)dz -» -oo (e-fQ). 
o 

We remark that even for a G C*(R) \Lip(R) the Garding inequality may fail as well: for example, 
with a(x) := -xlog \x\p{x) we have a G C£(R) n for all q G [1, oo[, but 

l l 

(Qv e , w £ )o = — 2 y p{ez)z log |ez|//(z)/9(z) dz < 2|loge / zp(ez)p' (z)p(z) dz —* — oo, 
o o 

since lim e ^o Jn z p( £z )p'( z )p( z ) dz = p(0) J„ zp'(z)p(z) dz < 0. 

Remark 2.2. (i) Let Q* denote the formal adjoint of Q with respect to the L 2 inner product (on 
x-space). Due to our regularity assumptions on a and c we have for any <p G H 1 (since a is real) 

TO 

Q*ip = *^2(-ajd Xj ip) + (c - div x {a))ip, 
i=i 

where the new coefficients —a, respectively c — div x (a), in place of a, respectively c, satisfy the 
exact same regularity assumptions, including the condition 



±div x (-a)-(c-div x (a)) = ^d?±_ c eL\[0,T};L^(M. n )). 

Thus the basic energy estimate (|2.16j) applies to ±9 t + Q* as well. In particular, the function h 
in the exponential factor occurring in the energy estimates is the same for Q and Q*. 

(ii) Although the method of derivation discussed above relied on a Garding-type inequality, it 
seems that in essence energy estimates are, in a vague sense, a necessary condition for a hyperbolic 
equation to hold in any meaningful context of "suitable Banach spaces of distributions" . In other 
words, whenever a hyperbolic differential equation can be interpreted directly in terms of such 
Banach spaces it allows to draw consequences on combinations of corresponding norms of any 
solution. For example, if the operator Q above generates a strongly continuous evolution system 
on some Banach space, then basic norm estimates for solutions follow from general principles of 
that theory (cf. [251155]). 

On the other hand, energy estimates are widely used to establish existence of solutions to (12.131) 
by duality and an application of the Hahn-Banach theorem. We recall the basic steps of such 
method in the following. 

2.1.2 Existence proof based on the energy estimate 

Let Rt := {(t,x) G R" +1 | t < T}. By abuse of notation we denote the trivial extension of a 
function v G C™{R T ) by zero for t > T again by v. Then C := {/ G C°°([0,T] x MP) \ 3v G 
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C™(R T ) with / = {-d t v + Q*v) | [0 ,t]xR«}- For < t < T and v G C c °°(i? T ) we use the notation 
w(t) := v(T - t) and g(t) := (-d t v '+ Q*v)(t). Then we have 

(d t +Q*(T-t))w(t) = g(T-t) 
w{0) = 

and an application of (|2.16|) (with Q* in place of Q; cf. Remark l2.2f i) above) yields 



sup \\w(r)\\ < 2exp(/ h(a)da) / \\(-d t + Q*v)(T - r)\\ dr = C h \ \\g(r)\\ dr. 

re[0.T] J J 



We may deduce that for / G £ x ([0, T];L 2 (M. n )) and v G C^{R T ) 

T T 

(f(r),v(r)) dr + {u , v(0)) < [ ||/(r)||o||«(r)||o dr + \\u \\ \\v(0)\\ 



<C sup \\w(r)\\ <CC h / ||ff(r)||odr, 
re[a,T] J 

where C depends on / and uq . Therefore the assignment g = (—dtv+Q*v) |[o,T]xR' it— v fo{f{ r )T v ( r )}odr+ 
{uq, v(0))o defines a conjugate-linear functional v: £^Con the subspace C of L 1 ([0, T]\ L 2 (R n )) 
such that < supo< r <T ||<7(?-)||o- Hahn-Banach extension of v yields a conjugate-linear func- 

tional v' : L 1 ([0, T]; L 2 (M. n )) — > C with the same norm estimate. 

Since L 1 ([0,T];L 2 (R n ))' = L°°([0, T];L 2 (R n )) there is u G L°°([0,T}; L 2 (R n )) such that u'(g) = 
(u,g) for all g G ^([0, T]; L 2 (tt")). When applied to 5 = (-ft<u + Q*u) | [0 ,r]xE» with w G C c °°(i? T ) 
we obtain 

T 

(2.17) J (u(t),-d t v(t) + (Q*v)(t)) dt=(u,((~d t + Q*)v) | [0 ,t]xM") 



T 

(f(t),v(t,.)) dt+ {u o ,v(0)) o . 



2.1.3 Model discussion of the weak solution concept 

Case of smooth symbol: If the coefficients of Q (and thus of Q*) are C°° then the above 
identity implies that u is a distributional solution to the partial differential equation Pu = / in 
]0, T[xR n . In fact, with <^ G C c °°(]0, T[xR n ) in place of v we have 

r 

((ft + Q)u, v ) = («, (-ft + Q*)(/3) = / (f(t), cf>(t, .)> di = </, p). 



Moreover, since Qu G L°°([0, T};H 1 (R™)) the differential equation implies that 

d t u = f-Qu 6i 1 ([0,T|;H- 1 (R B )) 

and thus u G AC([0, T]; ii" _1 (R n )). In particular, it makes sense to speak of the initial value 
u(0) G D'(R n ). Integrating by parts on the left-hand side of (|2.17[) (now reading (|2. 17[) from right 
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to left, and duality brackets in appropriate dual pairs of spaces) yields for any v G C^°(Rt) 

T T 

J (/(i),«(*))o dt + (tio, w(0)) = / (d t u(t) + Qu(t),v(t)\ dt - HT),v(T)) + («(0), v(0)) , 

</(t),f(t,.))o =0 

hence it(0) = uq 

Of course, uniqueness of the solution as well as more precise regularity properties can be deduced 
in case of C°° coefficients: For any s G R, / G L 1 ([0, T]; H s (R n )), and u G H s (« n ) the solution 
u is unique in the space C([0, T]; H s (R n )) (cf. [HI Theorem 23.1.2]). 

Case of non-smooth symbol: The weaker regularity assumptions made above imply Q*v G 
L 1 ([0,T];£ 2 (R n )) for all u G C™{R T ). We may thus define Qu £ X>'(]0, T[xR") by putting 

(Q« j¥ >) := (u\Q*<p)a V<p G C c °°(]0,T[xM"). 

Then equation (|2.17|) can be read as an equation in £>'(]0, T[xR n ), namely 

+ Qti, if) - (/, v>> G C c °°(]0, T[xR"). 

Furthermore, we can again show that the inital datum is attained: Note that in Qu — ^2 ajd Xj u 
each term can be interpreted as a multiplication of functions in L 1 ([0,T];H 1 (M. n )) with distribu- 
tions in L°°([0,T];H- 1 (W 1 )) (since u G L°°([0, T]; L 2 (R n ))) in the sense of the duality method 
(cf. [HI Chapter II, Section 5]). Applying Proposition 5.2 in 28J to the spatial variables in the 
products then yields Qu G L 1 ([0,T];W~ 1 ' 1 (M. n )). Reasoning similarly as above, the differential 
equation then gives 

d t u = f-Qu Gl^O.T];^- 1 - 1 ^")), 
which implies u G AC([0, T] ; W -1 ' 1 (K")) and further also that u(0) = u . 

Again higher regularity of u with respect to the time variable, namely u G C([0, T]; L 2 (R™)) can 
be shown by means of regularization and passage to the limit (e.g., similarly as in [2j proof of 
Theorem 2.8]). 



2.2 Regularization and energy estimates 

Several advanced theories make use of regularization techniques or concepts at crucial steps in 
their construction of solutions. Some of these theories succeed by regularization and a careful 
passage to the limit via energy estimates (as with Hurd-Sattinger and Di Perna-Lions theories 
presented below) . Others even base their solution concept on a further generalization of the weak 
solution concept beyond distribution and measure spaces and still obtain existence of solutions 
essentially from asymptotic stability of energy estimates (cf. the Lafon-Oberguggenberger theory 
below) . 

We introduce the following notation for partial differential operators that will be used in the sequel 

n 

(2.18) Pu := d t u + Ofc d Xk u + cu 

fe=i 

n 

(2.19) Lu := d t u + d X] (a,j u) + bu. 

3=1 
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2.2.1 Hurd-Sattinger theory 



We give a brief summary of the results from the first part in Hurd-Sattinger's classic paper [23:. 
We consider the Cauchy problem for the operator L on the closure of the domain Q, :=]0, oo[xR™. 

Definition 2.3. Let f 6 Lf oc (fl) and a,j (j = 1, . . . ,n) as well as b £ L; oc (fi). A weak solution in 
the sense of Hurd-Sattinger of the partial differential equation 

Lu = f on 

with initial condition uq £ L\ oc is a function u £ Lf oc (Cl) such that for all (j) £ C\ (K™ +1 ) we have 

/n 
( — u(t, x)(j)(t, x) — ^ ctj (t, x)u(t, x)d Xj 4>(t, x) + b(t, x)u(t, x)(f>(t, x)j d(t, x) 
- J =1 

f(t,x)cf)(t,x)d(t 7 x)+ / uq(x)4>(0, x) dx. 



n «*" 

Note that if all coefficients are C°° functions then a solution in the above sense solves the partial 
differential equation on in the sense of distributions. 

Theorem 2.4. Let aj (j — l,...,n), b, and f belong to Lf oc (fl) and uq £ Lf oc . Assume, in 
addition, that the following conditions are satisfied: 

(i) There exists C\ > such that for almost all (t, x) £ Q: ak(t, x) < C\ (k = 1 . . . , n). 

(ii) There exists a function /i £ L\ oc (\Q, oc[) ; /i > 0, such that b(t,x) > —^(t) for almost all 
(t,x) 6 0. 

(Hi) For each k £ {1, . . . , n} there exists < /ife G ^io C ([0, oo[) such that for almost all (t, x) £ fl 
a k (t,x)-a k (t,xi,...,x k -i,r,x k+ i...,x n ) 

> —f^k\t) f or almost all r £ JR. 

Xk - r 

Then there exists a weak solution u £ Lf oc (fl) to Lu — f with initial condition uo. 
Concerning the meaning of condition (iii) in Theorem 12.41 we mention two aspects: 

• In one space dimension we obtain gvgiLh^gfclJ > —^(f^ which resembles a one-sided Lip- 
schitz continuity condition in the x variable (apart from the fact that need not be 
finite or defined for all t). In particular, it excludes jumps downward (seen when going from 
smaller to larger values in the x argument). 

• Heuristically — replacing difference quotients by partial derivatives — condition (iii) can 
be read as diva(£,cc) > — ^2^k(t), thus giving an L 1 lower bound on the divergence of a. 
We observe that upon formally applying the Leibniz rule in the operator L we cast it in the 
form P as in ()2. 18(1 with c = div a + b. In combination with condition (ii) of Theorem 12.41 
we obtain that ^diva — c = — (^diva + b) has an L 1 upper bound (uniformly in x), which 
can be considered a substitute for condition (|2.14p used in the derivation of direct energy 
estimates in Subsection 2.1. 

Remark 2.5. Hurd-Sattinger ([H]) also give a uniqueness result for first-order systems in case 
of a single space variable and 6 = 0. For scalar equations the hypotheses require condition (i) to 
be strengthened to boundedness from above and from below and condition (iii) to be replaced by 
a Lipschitz property with an upper bound instead; in particular, no jumps upward are possible. 

Example 2.6. For the operator L in one space dimension and coefficients a(x) = sign(x) and 
b = 0, the Poupaud-Rascle theory is not applicable (as mentioned in [5D1 Section 1, Example 2]), 
but Hurd-Sattinger theory ensures existence of weak solutions, if the initial value belongs to £f oc . 
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2.2.2 Di Perna-Lions theory 



The weak solution concept introduced by Di Perna-Lions in [10] for the Cauchy problem for the 
operator P on a finite-time domain [0, T] x R™ can be interpreted in the following way. 

Definition 2.7. LetT > 0, 1 < p < oo, ±+± = 1, / G L 1 ([0,T\;L p (R n )), a k G L 1 {%T\-, L\ Q jR n )) 
(k = l,...,n), andce L 1 ([0, T] ; i^ oc (R n )) such that 

div(o)-c ei'dO.T];^^)). 

^4 function u G L°°([0, T]; L P (R™)) is called a weak solution in the sense of Di Perna-Lions of the 
partial differential equation 

Pu = f on]0,T[xR n 

with initial value Uq G L p (W l ), if 



(2.21) 




holds for all <p G C°° ( [0, T] , K n ) wit/i compact support in [0,T[xR n . 

Clearly, in case of C°° coefficients we obtain a distributional solution of the partial differential 
equation in ]0,T[xM™. 

Theorem 2.8. Existence of a weak solution u G L°°([0, T]; L p (IR n )) m f/ie sense of and with 
assumptions as in Definition \2. 7| is guaranteed under the additional hypothesis 

-div(a)-ce L 1 ([0,T};L oc (R n )), if p > 1, 
P 

div(a),c Gi 1 ([0,T];i oo (R")), if p = 1. 

Remark 2.9. Uniqueness holds in general under the additional hypotheses that c, div(a) G 
^([0, T];L°°(R n )), and for j = 1, . . . ,n also ^ G ^([0, T]; ^' C 9 (R™)) as well as 

G L 1 ([0,r]; J L 1 (M")) + L 1 ([0,T];L oc (M")). 

Example 2.10 (Hurd-Sattingcr applicable, but not Di Perna-Lions). Note that with a single 
spatial variable boundedness of div (a) = a' implies Lipschitz continuity. Hence, if a £ H (M) 
is not Lipschitz continuous but satisfies the one-sided Lipschitz condition in Hurd-Sattinger's 
existence Theorem 12.41 (iii), then a weak solution in the sense of Hurd-Sattinger to the problem 

d t u + d x (au) = / G L 2 (M 2 ), u | t=0 = G L 2 (R) 

is guaranteed to exist, whereas the general statement of DiPerna-Lions' existence theory (Theorem 
12.81 with p = q = 2) is not applicable to the formally equivalent problem 

d t u + ad x u + a'u = f G L 2 (R 2 ), u | t=0 = u G L 2 (R). 
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Example 2.11 (Di Perna-Lions applicable, but not Hurd-Sattinger). Let < a < 1 and consider 
the identical coefficient functions a± = a 2 G CT comp (R 2 ) (i.e., compactly suported functions in 
C^(R 2 )) given by 

ai(x,y) = a 2 (x,y) = -~(x - y)+x( x ,v), 

where x G T3(M. 2 ) such that x = \ near (0,0). Note that a\ is not Lipschitz continuous, since for 
x > but x sufficiently small the difference quotient 

ai(a:,0) -ai(0,0) _ x*- 1 
x a 

is unbounded as x — > 0. In particular, the latter observation shows that the Hurd-Sattinger 
existence theory is not applicable (condition (iii) in Theorem l2.4l is violated) to the Cauchy problem 
for the operator 

Lu = dt,u + d x (aiu) + d y (a,2u). 

On the other hand, we can show that with a = (01,02) the DiPerna-Lions existence theory is 
applicable to the Cauchy problem 

d t u + atd x u + a 2 d y u + (diva) u = f G L x ([0,T];i p (R 2 )), u | t=0 = u G i p (R 2 ). 



To begin with, we observe that 



x{x,y) 1 



dxai{x,y) = d x a 2 (x,y) = '-^ (x - y)+d x x(x,y) 

d y a x {x,y) = d y a 2 (x,y) = X X ^ 1 _ a (x - y)+d y x{x,y) 



div a(x,y) = --(x-y)1dW X (x,y) G C, a comp (R 2 ). 



yields 

1 
a 

Note that in the notation of Definition 12.71 and Theorem 12.81 we have c = diva G L°°(R 2 ) (and 
time- independent). Therefore, the basic assumptions for the solution concept to make sense as 
well as the hypotheses of the existence statement are clearly satisfied. 



As for uniqueness, we remark that all the conditions mentioned in Rcmark l2.9l are met if and only 
if a > 1/p. 

Remark 2.12. We mention that with coefficients as in the above example, the system of char- 
acteristic differential equations has forward-unique solutions, hence the Poupaud-Rascle solution 
concept for measures is also applicable. 



2.2.3 Lafon-Oberguggenberger theory 

The theory for symmetric hyperbolic systems presented in [25] by Lafon-Oberguggenberger allows 
for Colombeau generalized functions as coefficients as well as inital data and right-hand side. Thus 
we consider the following hyperbolic Cauchy problem in R n+1 

n 

(2.22) Pu = dt.u + a,jd Xj u + cu = f 

4=1 

(2.23) u | t= o= u Q , 

where aj (j = 1, ... ,71), c are real valued generalized functions in Q(M. n+1 ) (in the sense that all 
representatives are real valued smooth functions), / G Q(M. n+1 ), and initial value uq G CJ(R"). 
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The coefficients will be subject to some restriction on the allowed divergence in terms of e- 
dependence. A Colombeau function v G Q(R d ) is said to be of logarithmic type if it has a repre- 
sentative (v £ ) with the following property: there are constants N G N, C > 0, and 1 > r) > such 
that 

C 

sup \v e {y)\ < iVlog (— ) < e < r\ . 

j/GR d £ 

(This property then holds for any representative.) By a suitable modification of [271 Proposition 
1.5] it is always possible to model any finite order distribution as coefficient with such properties 
(in the sense that the Colombeau coefficient is associated to the original distributional coefficient) . 

Theorem 2.13. Assume that a,j and c are constant for large \x\ and that d Xk dj (k = 1, . . . ,n) 
as well as c are of logarithmic type. Then given initial data uq G G(W n ) and right-hand side 
f G g{R n+1 ), the Cauchy problem fU^ - fUg] has a unique solution u G g(R n+1 ). 

We also mention the following consistency result which shows that Colombeau theory includes 
the classically solvable cases: If we assume that the coefficients a,j and c are C°° then we have the 
following consistency with classical and distributional solutions (cf. [25] ) 

• If / and uq are C°° functions then the generalized solution u G C*(R n+1 ) is equal (in Q) to 
the classical smooth solution. 

• If / G L 2 (R;H s (R n )) and u Q G H s (R n ) for some s G K, then the generalized solution 
u G C/(R™ +1 ) is associated to the classical solution belonging to C(R; H S {W 1 )). 

Example 2.14. Consider the (1 + l)-dimensional operator 

Lu = d t u + d x (H(—x)u). 

Since the coefficient (of the formal principal part) has a jump downward neither Hurd-Sattinger 
nor Di Perna-Lions theory is applicable. In fact, it has been shown in [5D] Section 2] that none 
of the distributional products from the coherent hierarchy (cf. [28j and the introductary section) 
applied to H(—x) -u is capable of allowing for distributional solutions of the homogeneous Cauchy 
problem for arbitrary smooth initial data. 

Recall from Section 1 that measure solutions according to Bouchut- James exist for the corre- 
sponding Cauchy problem, if the Heaviside function (usually understood as a class of functions 
in L°°) is replaced by the particular Borel measurable representative with value at x = 0. For 
example, the initial value uq — 1 then yields the measure solution u = 1 + tS(x) in the sense of 
Bouchut- James as seen in Example 11.171 

However, Colombeau generalized solutions are easily obtained — even for arbitrary generalized 
initial data — if the coefficient H(—x) is regularized by convolution with a delta net of the form 
p £ (x) = \og{l/e)p{x\og{l/e)) (0 < e < 1), where p G C C °°(R) with j p = 1. Let a denote the 
class of this regularization in the Colombeau algebra Q, then the operator L may now be written 
equivalently in the form 

Pu = dtu + ad x u + a u, 

where a' w S and u G Q. Due to the logarithmic scale in the regularization the hypotheses of 
Theorem 12 . 131 are satisfied and the corresponding Cauchy problem is uniquely solvable. Moreover, 
for most interesting initial data (e.g. Dirac measures or L\ oc ) weak limits of the Colombeau solution 
u are known to exist and can be computed (cf. [20| Section 6]). In particular, for the initial value 
Uq = 1 we obtain the measure solution u = 1 + tS(x) as such a distributional shadow. 

Remark 2.15. (i) The basic results of Lafon-Oberguggenberger have been extended to the case 
of (scalar) pseudodifferential equations with generalized symbols in (TB]. Special cases and very 
instructive examples can be found in |26] , and an application of Colombeau theory to the linear 
acoustics system is presented in [27 J. 
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(ii) Colombeau-theoretic approaches allow for a further flexibility even in interpreting distribu- 
tional differential equations with smooth coefficients. For example, in [6] the concept of regularized 
derivatives is used, where partial differentiation is replaced by convolution with the correspond- 
ing derivative of a delta sequence. When acting on distributions this concept produces the usual 
differential operator actions in the limit. When considered as operators in Colombeau spaces, one 
can prove (cf. Theorem 4.1]) that evolution equations with smooth coefficients all whose deriva- 
tives are bounded have unique generalized function solutions for initial data and right-hand side 
in generalized functions. In particular, famous examples like the Lewy equation become solvable 
and Zuily's non-uniqueness examples become uniquely solvable then. 

2.3 Paradifferential techniques 
2.3.1 Energy estimates 

Bony's paradifferential calculus has been successfully applied in nonlinear analysis and, in par- 
ticular, to regularity theory for nonlinear partial differential equations. An ingredient in such 
approaches is often a refined regularity assessment of corresponding linearizations of the differen- 
tial operators involved. A recent account of Metivier's methods and results of this type can be 
found in Subsection 2.1.3], or with more details on microlocal properties in [17] . 

Let s£K and H^(R n ) denote the Sobolev space H s (R n ) equipped with the weak topology. We 
consider a differential operator of the form 

n 

P v (x,t;d t ,d x ) := d t + y^aj{v(x, t))dj, 
j=i 

where a 3 G C°°(R) (j = l...,n) and v G L°°([0, T]; H s (R n )) n C([0, T]; H^(R n )) such that 
d t v e L°°([0,T];fl' s - x OR n ))nC([0,T];ir»- 1 (R n )). 

Remark 2.16. Not all hyperbolic first-order differential operators with coefficients of regularity 
as above can be written in the special form of P v . In fact, this amounts to writing any given 
list wi, . . . , w n of such functions as wj = aj o v (j = 1 . . . , n) with a,j G (7°°(R) and v as above. 
The latter is, in general, not possible, which can be seen from the following example: consider the 
Lipschitz continuous functions wi(t) — \t\ and w 2 (t) = t; if W\ = a\ o v and w 2 = a,2 ° v with a 
Lipschitz continuous function v, then v is necessarily non-differentiable at 0; on the other hand 

1 = «4(0) = lim(a 2 (v{h)) - a 2 (v(0)))/h = lim a' 2 (£(h))(v(h) - v(0))/h, 

where £(h) lies between v(0) and v(h); hence a' 2 (^(h)) — > a' 2 (w(0)) and the second factor (v(h) — 
v(0))/h stays bounded, but is not convergent; in case a' 2 (v(0)) — we obtain the contradiction 
1 = 0, in case ^(^(O)) ^ we have a contradiction to convergence of the difference quotient for 
w 2 . 

The key technique in analyzing the operator P v is to replace all terms dj(v)dj by T a .r v )dj, i.e., 
partial differentiation followed by the para-product operator T a .r v ), and then employ estimates 
of the error terms as well as a paradifferential variant of Garding's inequality (cf. 2, Appendix 
C.3-4]). This leads to the following result. 

Theorem 2.17 ( 2, Theorem 2.7]). If s > | + 1, then for any f G L°°({0, T]; iP(R n )) n 
C{[0,T};H^(R. n )) and u Q G H s (R n ) the Cauchy problem 

P V U = /, U \ t = 0= UQ 
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has a unique solution u G i 2 ([0, T]; H s (R n )). Moreover, u belongs to C([0, T]; H s (R n )) and there 
are constants K, 7, C > suc/i i/iai u satisfies the energy estimate 

t 

\\u(t)\\ 2 s <Ke^\\u(0)\\ 2 s +C J e^-^WPMr^dr. 



2.3.2 Improvement of regularity in one-way wave equations 

We briefly recall some basic notions and properties concerning symbols with certain Holder regu- 
larity in x and smoothness in £ a la Taylor (cf. [35]). 

Definition 2.18. Let r > 0, < 5 < 1, and m G R. A continuous function p : R n X R n — * C 
belongs to the symbol space C£S , j Tl (5 , if for every fixed x G R™ £/ie map £ 1— » p(x,£) is smooth and 
for all a G Nq there exists C a > such that 

\d? P (x,0\ < C a (l + |£|) m ~ |Q| Vx,£ g R™ 

and 

II^pCOIIcj < c a (i + leir-i"^ v£ e R n . 

Basic examples are, of course, provided by symbols of differential operators »a<9 Q with coefficient 
functions a a G C£ (|a| < m) or any symbol of the form p(x,£) = a(a;)/i(a;, £), where a G C£ and h 
is a smooth symbol of order m. 

Symbol smoothing: By a coupling of a Littlewood-Paley decomposition in £-space with convo- 
lution regularization in s-space via a (5-dependent scale one obtains a decomposition of any symbol 
p G CIS™ S in the form 

p = p 8 + p\ where p s G Sft and p b G CIS^ 5 . 

Observe that p" is C°° and of the same order whereas p b has the same regularity as p but is of 
lower order. 

Mapping properties: Let < 5 < 1 and —(1 — S)r < s < r. Then any symbol p G C^S™ S 
defines a continuous linear operator p(x,D) : H s+m (M. n ) — > £P(R"). 

Elliptic symbols: p G CIS™$ is said to be elliptic, if there are constants C, i? > such that 

\p( X ,o\ >c(i + ieir veGK",i^i>i?. 

One-way wave equations result typically from second-order partial differential equations by a 
pseudodifferential decoupling into two first-order equations (cf. [SH Section IX. 1]). For example, 
this has become a standard technique in mathematical geophysics for the decoupling of modes in 
seismic wave propagation (cf. [32] ) . The corresponding Cauchy problem with seismic source term 
/ G C oo ([0,T];H s (W 1 )) (with s G M) and initial value of the displacement u a G H s+1 {W l ) is of 
the form 

(2.24) d t u + iQ(x,D)u = f 

(2.25) u | t _ = u , 

where Q has real-valued elliptic symbol q G C^S 1 with r > s. 
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Lemma 2.19. If q 6 C r S™ is elliptic, then q$ £ S™ s is also elliptic. 

Proof. By ellipticity of q and the symbol properties of q b there are constants C\, C%, R > such 
that 

Cx(i + iei) m < |g(*,OI < \qH^O\ + < + + l£ir~ r<5 

holds for all with |f| > i? > 0. Therefore 

>(Ci-^(i + |^|)- r5 )(i + ieir>C(i + |e|) m Vx^eMMfl ># 

for suitably chosen constants C and R' > 0. □ 

Let < <5 < 1. We have the decomposition g = g" + 9 b , where g" € 5j a and q b E C r S\^ s Sr . By 
Lemma \2. 191 = g"(x, D) is elliptic and thus possesses a parametrix £ S^. 

We have 

(flt + iQ)E*f = {d t + iQ 1 + %Q v )E^f = d t E*f + iQ*E*f + iQ b E*f = d t E*f + f + iR?f + iQ^E^f, 
where R* is a regularizing operator. Therefore 

(fit + iQ){u - E*f) = -dtEflf - iR*f - iQ b E*f =: /, 
where the regularity of the right-hand side / can be deduced from the following facts 
d t E^f G C°°([0,r];ff s+1 (R n )), R^f £ C°°([0, T}; H°°(W 1 )), Q b E t f 6 C*°°([0, T]; H s+Sr (W 1 )). 
Hence / € C*°°([0, T]; H s+min ( Sr M(R n )). 

If we put w = u + E$ f and wo := + E^f(0), then the original Cauchy problem (|2.24II2.25|) is 
reduced to solving the Cauchy problem 

d t w + iQ{x,D)w = f, w\ t =o=wo, 

where the spatial regularity of the source term on the right-hand side has been raised by min((5r, 1). 

Remark 2.20. In case of a homogeneous (1 + l)-dimensional partial differential equation the pre- 
cise Holder-regularity properties of classical as well as generalized solutions have been determined 
in QU Section 3]. 
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